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CORRELATION OF FATIGUE CRACK GROWTH RATE AT DIFFERENT 
STRESS RATIOS FOR QUENCHED AND TEMPERED STEELS AND 

OTHER ALLOYS 

I. M. Robertson 

Department of Defence. DSTO Materials Research Laboratory. PO Box 50, Ascot Vale. Vic. 3052. Australia 

Abstract —Measurements of the effect of stress ratio on the constant amplitude fatigue crack growth rates 
in four quenched and tempered steels in the Paris regime arc reported. This data and published data for 
other alloys (including lower strength steels and non-ferrous alloys) arc evaluated, and a correlation 
function suitable for practical fatigue life calculations is derived. In addition to stress intensity factor 
range and stress ratio, other significant parameters arc the yield stress of the material and its thickness. 

For the four steels on which new measurements were made, the degree of dependence of the crack growth 
rate on stress ratio may be related to sensitivity to environmental conditions. 


NOMENCLATURE 

a = crack length 
C. D = Paris coefficients 
AK = full stress intensity factor range 
AK,, = effective part of A A 
A K, m = tensile part of A K 
n = Paris exponent 
.V = number of fatigue cycles 
= maximum load in fatigue cycle 
R = stress ratio 
i = thickness 

(/(£)- crack closure function 
P(RI = crack growth rate correlation function 
i7y yield stress 


INTRODUCTION 

The present work is concerned with fatigue crack growth rates in the Paris regime (or Stage II 
crack growth region). To a first approximation, fatigue crack growth in the Paris regime in metallic 
materials under constant amplitude loading is described by the Paris-Erdogan relation: 

da 

— = OAK)" (1) 

d N 

There have been a large number of studies resulting in the modification of Eq. (1) to incorporate 
the second order effect of mean stress (or stress ratio) and other variables. Some of the equations 
that have been proposed were tabulated by Chand and Garg [1]. 

In the present work the formulation presented by Elber [2] is followed. The crack growth rate 

327 


f 







.128 


I M RoBHtrsoN 


is regarded as a function of AK and the stress ratio. R: 


du 
d N 


du du du 


12 ) 


where U is regarded as a function of R only. Although Elber associated U{R) with the portion of 
AK for which the crack is actually open, in the present work the function is treated simply as an 
empirical function for the correlation of fatigue crack growth rates in a material under different 
stress ratios. As discussed by Schijve [3], this type of lHR) can be obtained from experimental 
data by measuring the relative displacement of logfdu/dAf) versus log(AK) curves for different 
stress ratios. 

The Paris-Erdogan relation is rewritten in the following way to incorporate stress ratio: 
du 

— = D(AK efr ) n = D[l/(R)AK] B = C(R)(AK)" (3) 


The growth rates under different stress ratios at the same AK arc then: 


^ (AK,. R,) = D[t/(R, )AK, ]" 
diV 


= [C/(R, )/(/{«„)]" ~(AK,.R„) (4) 

For equal growth rates at different stress ratios the stress intensity factor ranges required are given 
by: 


AK 0 = AK, [!/(/?, )/(/(/?„)] (5) 

It should be noted that these equations apply only within the Paris regime and cannot be used 
for the threshold region or the high crack growth rate region (Stages I and III), where the effect 
of stress ratio is more pronounced. 

Several different forms have been proposed for the function U(R). Some of them arc as follows: 

(a) A basic correction which converts the full stress intensity factor range. AK. into the tensile 
part of the range, AK lent . This will be referred to as the "basic correction", and is commonly the 
only correction for stress ratio effects that is applied in fatigue life calculations for structures 
containing residual stresses and other sources of mean stress. 


U = I for R > 0 

= 1/(1- R) for R < 0 


( 6 ) 


(b) Elber's formula (2): 


(/ = 0.5 + 0.4R for — 0.1 < /? <0.7 (7) 

(c) Schijve’s formula (Finney and Deirmendjian [4]) for aluminium alloys in the range 
-1 <R<0.7 


U = (0.55 - 0.2R - 0.25R 2 - 0.1 R 3 )/( 1 - R) 


( 8 ) 


and similar polynomial expressions given in Refs [3] and [5]. 
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(d) The formula of Kurihara t*r at. [6] for pressure vessel steels: 

U = 1 for 0.5 < R < O X 

(9) 

= I(1.5-X) for -5 <R <0.5 


(e) The formula of Eason et at. [7] for ferritic steels in the range 0 < R < 0.9 

V = 1 /(2.88 - X) (10) 

If the da/dN versus AX curve is known for a single stress ratio. X 0 , for the material in question 
(usually /? o = 0), the steady state growth rate under any loading AX,. X, can be calculated from 
Eq. (4). This of course requires knowledge of )/C/(f? 0 ). Conversely, the value of U(R,)/U(R 0 ) 
can be obtained from fatigue crack growth rate curves for R = R, and R = R„. It is not possible 
to obtain absolute values of U(R) in this way, but neither is it necessary to have absolute values 
of U(R ) in order to calculate growth rates for stress ratios other than those for which data are 
available. AH that is necessary is the ratio t/(X, )/l/(X 0 ). 

We could use the function: 


V(R) = U(R)/V(0) (11) 

as an empirical fatigue crack growth rate correlation function. However, it is more instructive to 
incorporate the basic correction (Eq.(6)) from the outset. This results in a correlation function 
identified here as K'(K). It is derived as follows. 

Equations (1)—(5) are modified by replacing AX with AX len , (including new values of the 
material constants C and D) and the functions U{R) and C{R) are replaced respectively by the 
functions: 


l/’(X)=(/(X) 

for X > 0 

(12) 

= (l -R)U(R) 

for X < 0 

cm = cm 

for X > 0 

(13) 

= (I-X)C(X) 

for X < 0 



The empirical fatigue crack growth rate correlation function to be used here is then: 

nx) - ummo) 


= V(R) for XX) 


(14) 


= (1 — X)K(X) for X < 0 

From Eq. (5), the crack growth rate at stress ratio X would equal that at R = 0 when the tensile 
component of the stress intensity factor range is: 

AX leo ,. R = AX, en , 0/ r(X) (15) 

where AX tefl ,, 0 is the tensile component of the SIF range at X = 0 (equal to AX). From Eq. (4) the 
growth rates at constant AX tent are: 

da da 

— <AX, en ,. R ) = — (A x ten , 0 )[nx)] n (16) 

Figure I shows the function K'(X) derived from Eqs (7)-( 10). (12) and (14). Although the U{R) 
functions differ quite significantly, the values of the corresponding F'(X) functions are similar 
because of the normalisation to the value U{ 0). 
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Stress ratio 

Fig. I. Values of the crack growth rate correlation function l'(RI derived from tqs 1 7 > tot 10) for the 

closure function l IR). 


One of the goals of the present work was to decide why different materials or testing conditions 
give rise to different V'{R) correlation functions. Four quenched and tempered steels were tested 
over a range of stress ratios between - 2 and + 0.7 to examine the effect of yield strength on the 
correlation function F'(R). Published fatigue crack growth rate curves for other alloys were also 
re-examined and evaluated in terms of V'(R). 

MATERIALS AND METHODS 

The compositions and tensile properties of the steel plates from which specimens were machined 
are listed in Table 1. All are high toughness, quenched and tempered steels used for submarine 
hull and other applications. The strength depends on the composition and the temperature of the 
tempering treatment applied to the steels after quenching. Centre-cracked fatigue specimens with 
a width. H< of 100 mm and thickness. B, of 12.5 mm were machined, and stress relieved at 600 C 
before grinding to final dimensions. Fatigue testing was carried out in laboratory air (about 20‘ C 
and 50% RH) under sinusoidal loading at a frequency of 3-5 Hz. 

For each specimen the stress ratio was increased in a stepwise manner (minimising transient 
effects as described by Robertson [8]) so that the cracks extended by about 4 mm for each step 
of constant amplitude loading but the growth rate remained approximately 10* 7 m cycle. The 
loads that were applied are listed in Table 2. A single specimen of each steel was tested except for 
Q1N. which was used to confirm the reproducibility of the result [8]. Fatigue crack growth was 


Table I. Composition (V.wt) and tensile properties of the four steels 


Steel 

%C 

%Si 

%Mn 

%P 

%S 

%Ni 

%Cr 

%Mo 

%Cu 

•4 V 

QIN 

0.16 

0.25 

0.31 

0.010 

0.008 

2.71 

1.42 

0.41 

0.10 

001 


Proof stress: 570 MPa 

UTS: 663 MPa 






Q2N 

0.15 

0.35 

0.31 

0.010 

0.003 

3.20 

I 53 

0.45 

0.08 

0.01 


Proof stress: 707 MPa 

UTS: 823 MPa 






BIS 

0.14 

0.27 

0.98 

0.010 

0.001 

1.30 

0.52 

0.37 

0.21 

0.02 


Proof stress: 700 MPa 

UTS: 770 MPa 






HYI30 

0.10 

0.32 

0.79 

0.006 

0.003 

4.90 

0.54 

0.49 

0.05 

0.07 


Proof stress: 958 MPa UTS: 1004 MPa 
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Table 2. Sequence of loads applied to specimens U' m „ m kN) 


QIN 

R 

-2.03 

-1.04 

-0.51 

001 

n 29 



0.60 




109.0 

108.0 

109.0 

108.5 

1250 



196.0 


Q2N 

R 



-0.48 

0.01 

O20 

0.41 


0.60 

0.71 





116.5 

110.0 

10’ 5 

110.0 


144.0 

148 5 

BIS 

R 



-0.49 

0.01 

0 20 

040 


0 61 

0 70 


P~. 



125.0 

119.0 

121.5 

127.0 


167.5 

166.5 

HY 

R 

-1.98 

—1.01 

-0.49 

0.03 

0 28 


0.51 


0.71 


P~. 

122.5 

108.5 

98.0 

945 

103.0 


127.0 


159.5 


monitored using a travelling microscope, and growth rates calculated using the seven point 
incremental polynomial method. 


RESULTS 


Quenched and tempered steels 

Results for HY130 from the present wrork are shown in Fig. 2. where the full range of stress 
intensity factor has been plotted (to avoid overlapping curves). Results for Q1N were published 
previously (Fig. 3 of Ref. [8]). The remaining steels. Q2N and BIS812EMA. gave results similar 
to Q1N and HY130 except that testing was restricted to stress ratios between -0.5 and +0.7. 

The values of the function K'(R)= U'(R),U'(0) for the different steels are compared in Fig. 3. 
Some additional data points extracted from the relative displacement of the crack growth rate 
curves of Jones [9] are also included (for Q1N and HYI30 in laboratory air at R = 0.5 and a 
growth rate of 10“ 7 m/cycle). The sensitivity of the crack growth rate to stress ratio appears to 
increase with the yield strength of the steel at positive R, but the opposite trend applies for 
negative R- 

The data points shown in Fig. 3 are approximated by the equations: 


QIN 

R > 0: 

V’(R)= 1 + 0.16R 

(17) 

BIS 

R >0: 

r<R)= 1 +0.34R 

(18) 

Q2N 

R> 0: 

F'(R) = 1 + 0.36 R 

(19) 


to* 

£ 

y 

>» 

u 

£ 

2 

■o 

•D 

10 ' 


10 20 30 50 70 

&K, MPa^n 

Fig. 2. Fatigue crack growth rate curves for HY130 steel in air for the stress ratios. R. of -2. - I. -0.5. 
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Fig. 3. Values of the crack growth rate correlation function l '(RI for the steels 01N. 0-N. BISXI3EMA 
and HY130 from the present work Hilled symbols) and from Jones (9] (open symbols at K — 0.51. 


HY130 

R> 0: 

1"(R)= 1 + 0.44R 

(20) 



f F(R) = 1.25(1 -R)/( 1.25- R) 

(21) 

All 

R <0: 

(22) 

lv"(R) = 1 + log( 1 - R)/2.7 


Additional values of VfR) derived from published crack growth rate curves [6.10-13] for 
quenched and tempered steels are shown in Fig. 4. There are two sets of results for Q1N steel 
(both calculated from the relative displacement of growth rate curves at 10~ 8 m/cycle). Although 
James and Knott [10] found very little influence of R on the growth rate of long, through cracks 
in Q1N for R between 0.2 and 0.7 (similar to results for Q1N in Fig. 3). Cowling and Knott [11] 
found much greater sensitivity to R for short, semi-elliptical cracks in QIN (even greater sensitivity 
than for HY130 in Fig. 3). As James and Knott, and Cowling and Knott did not report results 
for zero R, an extrapolation has been used to position the data point at R = 0.2 in each case. For 
example, the data of Cowling and Knott allow the values V( 0.35)/K(0.2), 1(0.5) 1(0.2) and 
l'(0.7)/K(0.2) to be determined, and extrapolated to obtain F(0)/K(0.2). The reciprocal of this is 
1"(0.2) and the other values plotted are given by F'(R)= F'(0.2)[K(R)/F(0.2)]. 



Stress ratio 

Fig. 4. Values of the correlation function F7/0 derived from published fatigue crack growth rate curves 
for quenched and tempered steels. Values were obtained at growth rates between It) * and 10 ' m cycle. 
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Fig. 5. As for Fig. 4 but for steels with yield stresses between 2tX) and -MM) M Pa. 


Other steels 

Except for the quenched and tempered steels considered above (with yield strengths ranging 
from about 550 to 1000 MPa), relevant data for steels appear to be available only for yield 
strengths between 200 and 400 MPa. The values of the correlation function l'(K) derived from 
published crack growth rate curves [5,6,14-19] for these steels are shown in Fig. 5. 

There is a considerable amount of scatter, as might be expected with data derived from a wide 
range of sources. However, trends are still apparent. For positive stress ratios. V'{R) increases most 
rapidly with increasing R for lower strength steels tested at lower thickness. This can be attributed 
to larger plastic zone size relative to the thickness of the specimen, cracking under more-nearly 
plane stress conditions, and a greater degree of slant fracture (see Schijve [3]. McEvily [20] for 
discussion of increased crack closure under plane stress conditions). This is the opposite trend to 
that observed for the quenched and tempered steels, where essentially plane strain conditions 
prevail and the fracture faces are nearly flat. 

For negative stress ratios there are fewer data but it appears that steels of higher strength are 
less sensitive to R (as for the quenched and tempered steels). 

Austenitic stainless steel, aluminium and titanium alloys 

Figure 6 shows values of V'(R) derived from published V(R) functions for aluminium and 
titanium alloys [3], and from crack growth rate curves for an aluminium alloy [21] and type 304 
stainless steel [15,22]. Most of the data are for positive R and show a strong dependence of crack 
growth rate on stress ratio for R > 0. As for the low strength steels considered above, this can be 
attributed to the large degree of crack closure and slant fracture due to the low thickness at which 
most of the alloys were tested. 


DISCUSSION 

Examination of Figs 3 and 5 suggested that yield strength and specimen thickness influence the 
values of the correlation function. Figure 7 shows the effect of yield stress for both groups of steels 
considered above. For most of the steels, F'(A) increases approximately linearly as R increases or 
decreases from zero. Therefore, stress ratios of +0.5 and -0.5 have been selected to represent 
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Stress ratio 

Fig. 6. As for Fig. 4 but for aluminium alloys [3.211. i\pc 304 stainless steel ( 15.221 and a titanium 

alloy | 31 ‘ 



Fig. 7. Values of F'(0.5) and I f -0.51 plotted against yield strength for the steels shown in Figs 3, 4 

and 5. 


positive and negative R respectively in Fig. 7. In some cases it was necessary to interpolate I linearly) 
between data points to obtain the V'(R) value at R = 0.5 or -0.5. 

There is a large amount of scatter, especially for the lower strength steels. Some of this can be 
attributed to the fact that A/C, en „ R and n y (or thickness) arc not sufficient to characterize the 
material and crack growth conditions. However, a targe part of the scatter probably arises from 
inaccuracies in the experimental determination of growth rates at different R. and in extracting 
the V'(R) values from published crack growth rate curves. For the quenched and tempered steels 
there are two anomalous data points. The high value of F'(0.5) for semi-elliptical cracks in Q1N 
steel [11] has been discussed above. The high value of K'(-0.5) for a steel with a yield stress of 
794 MPa [12] can be attributed to experimental scatter, as the corresponding F'(-l) value is 
lower (Fig. 4). 

For negative R the effect of yield strength on F'(R) is weak, but F'(R) is slightly higher for the 
lower strength steels. Specimen thickness does not appear to be a significant influence. A possible 
mechanism for the yield stress dependence is the crushing together of asperities on the fracture 
faces during the compressive part of the fatigue cycle. Linear regression relations for the data are 
as follows: 
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Quenched and tempered: 

Q and T (excluding [12]): 
Lower strength steels: 
Combined data (not [12]): 


F'(-0.5) = 1.15-0.11 X 10 J (rr y /MPa) 
(6 points, r = -0.411 

K'(-0.5)= 1.16-0.14 x 10'' </r y /MPa) 
(5 points, r = -0.82) 

K'(-0.5)= 1.59- 1.4 x 10 J (ff y /MPa) 
(4 points, r = - 1.00) 

r(-0.5)= 1.10-0.06 x 10‘ J (tr y /MPa) 
(9 points, r = -0.48) 


(23) 

124) 

(25) 

(26) 


For positive stress ratios, the quenched and tempered steels and the lower strength steels behave 
differently, but the effect of yield stress is much stronger than for negative R for both groups of 
steels. Thickness also appears to be significant, but only for the lower strength steels. Table 3 shows 
the results of linear regression analyses for F'(0.5) against rr y . t. a y t and (1000 MPa - a y ) 2 . The 
parameter o$t was selected because the ratio of specimen thickness to crack-tip plastic zone size 
is proportional to 

'/(Kimu/tfy) 2 - Vyt/Kmn (27) 

and Km, is already determined by AK and R. The anomalous data point for semi-elliptical cracks 
in Q1N [11] has been excluded. 

It is recognized that crack closure at positive stress intensity is more prevalent under plane 
stress conditions than under plane strain [3,17,19,20.23], and that closure is largely irrelevant 
during the initial stages of fatigue cracking of a structure. On this basis we would expect that 
K'(A) for positive R would decrease towards unity as the thickness or yield strength of the material 
increased. This is the trend observed for the lower strength steels in Fig. 7 and Table 3 (and suggests 
that crack closure is the mechanism for the dependence of crack growth rate on stress ratio in 
these steels). 

However, the quenched and tempered steels do not behave in this way (either in the general 
trend of V'(R) with yield stress, or in the specific case of semi-elliptical cracks in Q1N [11] where 
strong R dependence was observed). A possible explanation for the greater dependence of crack 
growth rate on stress ratio for the higher strength steels lies in their greater sensitivity to environ¬ 
ment (rather than yield strength per se ). The steels Q1N and HY130 were examined by Jones [9] 
in dry air, laboratory air and sea water. The steels HY80 and HY130 (HY80 is almost identical 
to Q1N) were examined by Kwun and Fine [24] in laboratory air and dry argon. These investi¬ 
gations showed that the two steels have similar crack growth rates in dry air or dry argon, but 


Table 3. Linear regression equations for HO.5) for steels. Yield stress. <r r . 
in MPa, thickness, i. in mm 


Steel 

group 

No. of 
points 


H0.5) 

Regression 

coefficient 

Q and T 

10 

1.23 

-0.99 x !0- ft (l000-<r y ) 2 

-0.87 



0.81 

-t- 0.46 x I0' 3 (<j y ) 

+0.83 



1.06 

+ 0.014 x IO-' , (ff 2 t) 

+0 51 



1.30 

-0.013(f) 

-0.46 

Lower a y 

8 

1.81 

— 1.7 x I0' 3 (<7 y ) 

-0.65 



1.36 

— 0.0S 1 X IO *(<T?f) 

-0.69 



1.38 

-0.013(f) 

-0.71 
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Fig. 8. Schematic diagram for the effect of yield strength and thickness on the crack growth rate 

correlation function V'(R). 

the growth rate increases much more for HY130 when the environment is changed to laboratory 
air than it does for Q1N (the same trend does not apply for the seawater environment). 

The greater sensitivity of HY130 to moisture in the air could be expected to result in greater 
sensitivity to stress ratio at R > 0, where higher stress ratio opens the crack more and makes the 
crack tip more accessible to the environment. 

Some results are not explained by this argument and are included here for completeness: 

(a) On changing the environment from dry air to seawater. Q1N shows a larger increase in 
crack growth rate than HY130, but the sensitivity to R of the crack growth rate in HY130 in 
seawater is still higher than that of Q1N [9]. The values of K'(0.5) for both Q1N and HY130 are 
similar in seawater, laboratory air and dry air. 

(b) Cowling and Knott [11] observed strong dependence on R of the growth rate of semi¬ 
elliptical cracks in Q1N, in spite of weak environmental sensitivity (and plane strain conditions). 
This may be related to the fact that the cracks were short. 

Figure 8 summarizes the present discussion in schematic form. For negative R, the correlation 
function K'(R) decreases slightly with increasing yield stress but is not apparently affected by 
thickness. For positive R, there are two conflicting trends with increasing yield stress. For low 
yield strengths (and thicknesses) V’(R) decreases with increasing yield stress (or thickness) as 
conditions of plane strain and flat fracture are approached. The change in K'(/?) can be attributed 
to a crack closure mechanism. For high yield strengths (and thicknesses) closure is not significant, 
but V'(R) increases with yield strength because there is a concomitant increase in sensitivity to the 
testing environment. The correlation between yield strength and environmental sensitivity is 
unlikely to persist for steels whose compositions or microstructures differ widely from those of the 
quenched and tempered steels of the present study. Some of the scatter for the lower strength 
steels may be due to variations in their environmental sensitivity. 

CONCLUSIONS 

In this paper, measurements of the effect of stress ratio on the fatigue crack growth rate in 
quenched and tempered steels with a range of yield strengths have been reported. The results have 
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been compared with published measurements for other steels, aluminium alloys, austenitic stainless 
steel and a titanium alloy. 

A crack growth rate correlation function F'(R) has been derived which enables an experimental 
crack growth rate curve for R =0 (or any available R) to be converted into a curve for any other 
R. The function F'(R) is derived from du/dlV versus A K curves, so it does not specifically relate 
to any particular mechanism. 

This function reveals more clearly than previously proposed correlation functions (such as the 
crack closure function U) the effects of yield stress and thickness, because it reduces the amount 
of scatter. 

The approach used here is to gradually refine the crack growth rate prediction. The major 
influence is of course the tensile part of AK. Of secondary importance is R. Further refinement is 
introduced by adjusting the R correction (i.e. the correlation function F'(R)) to take yield strength 
and thickness into account. 

For steels, yield stress has a slight effect on crack growth rate for negative R. For positive R. 
the effect is greater but the growth rate becomes less sensitive to R as plane strain conditions are 
approached (with increasing strength or thickness), with the exception that higher strength 
quenched and tempered steels show increasing sensitivity to R. The latter trend is possibly due to 
increasing sensitivity to environment, but the evidence for this is not conclusive. 

Aluminium alloys, austenitic stainless steel and titanium alloys have not been considered in the 
same detail, but show strong dependence of the crack growth rate on R. This is attributed to the 
conditions of plane stress under which most of these materials have been tested. 
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